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ABSTRACT 

The transit signature of exoplanets provides an avenue through which characterization of exoplan- 
etary properties may be undertaken, such as studies of mean density, structure, and atmospheric 
composition. The Transit Ephemeris Refinement and Monitoring Survey (TERMS) is a program to 
expand the catalog of transiting planets around bright host stars by refining the orbits of known plan- 
ets discovered with the radial velocity technique. Here we present results for the HD 38529 system. 
We determine fundamental properties of the host star through direct interferometric measurements of 
the radius and through spectroscopic analysis. We provide new radial velocity measurements that are 
used to improve the Keplerian solution for the two known planets, and we find no evidence for a previ- 
ously postulated third planet. We also present 12 years of precision robotic photometry of HD 38529 
that demonstrate the inner planet does not transit and the host star exhibits cyclic variations in 
seasonal mean brightness with a timescale of approximately six years. 

Subject headings: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (HD 38529) 



1. INTRODUCTION 

The study of transiting exoplanets has undergone a 
remarkable evol ution since the first transit detection 
in HP 209458b (jCharbonneau et alj [20001 : IHenrv et all 
2000). Projects such as the Hungarian Autom ated Tele- 
scope Netw ork (HATNet) (iBakos et al.ll2004f) . and Su- 
perWASP (|Pollacco et all I2006D are routinely detect- 
ing new transiting planets. Results from the NASA 
Kepler mission are breaking new ground by discover- 
ing planets that are smaller and at longer orbital pe- 
riods than those acce ssible from ground-based surveys 
(jBorucki et al.ll2l)Tl a.b'). However, the most important 
transiting planets for follow-up observations and char- 
acterization of atmospheres continue to be those dis- 
covered first by the radial velocity (RV) method, be- 
cause of the bias towards bright host stars. The Transit 
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Ephemeris Refinement and Monitoring Survey (TERMS) 
seeks to provide more of thes e opp ortunites through 
orbital refinement (|Kane et alj 12009ft , particularly for 
plan ets in eccentric orbits wit h higher transit probabili- 
ties (|Kane &: von Braunll2008f) . TERMS also allows the 
study of long-term variability of the host star such that 
it can be correlated with the stellar mag netic activity cy- 
cle in the context of being a planet-host (jPragomi r et al.1 
I2012T) . 

Multi-planet systems present particularly interesting 
cases since studying the dynamical interaction of the 
planets is an additional advantage to refining the or- 
bits. One such TERMS target is the planetary system 
orbiting HD 38529, whose planets have been discovered 
via precision RV measurements. The inne r planet is in 
a ~ 1 4 day orbit and was discovered by iFischer et al.1 
(120011) . T he second plan et, in a ~ 6 year orbit, was dis- 
covered bv IFischer et al.l (|2003[ ). The system was further 
investigated in th e cont ext o f multi-planet syste ms by 
iWittenmver et all (|2009f ) and lWright et al.1 ((2009), both 
of whom pr ovided new RV data and revised orbits for 
the planets. iBenedict et al.l (|2010l ) furthered the studies 
of the system by providing new RV data as well as Hub- 
ble Space Telescope astrometry. Their analysis indicated 
that there may be an additional planet located at an or- 
bital period of ~ 194 days, though this was not confirmed 
by their observations and they encouraged further study 
to help resolve the issue. 

In this paper we present an exhaustive analysis of both 
the host star and the planets in the HD 38529 system. 
We performed interferometric observations of the host 
star using the CHARA Array to determine a direct mea- 
surement of the radius, which we compare with that de- 
rived from an analysis of high-resolution spectra. We 
provide new RV data for the system, improving the orbits 
of the planets and greatly extending the time baseline. 
Our analysis of the combined data finds no evidence for a 
third planet located near 194 days. From these new data, 
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TABLE 1 

Stellar. Properties of HD 38529 



Value Value 

Parameter Spectroscopic Interferometric 

6» UD (mas) 0.593 ±0.016 

6> LD (mas) ■•■ 0.611 ±0.016 

Luminosity (Lq) ■■■ 5.777 ±0.186 

Radius R t (R Q ) 2.34 ±0.07 2.578 ± 0.080 

T eS (K) 5619 ± 44 5576 ± 74 

[Fe/H] 0.38 ±0.03 

vsini (km s" 1 ) . 3.20 ± 0.50 

logg 3.83 ±0.06 

Mass M* (M ) . 1.36 ±0.02 
Age (Gyr) 4.45 ± 0.23 



Reference 



this work (\ 
this work ( 
this work ( 
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Note. — 

we calculate an accurate transit ephemeris and predict 
the parameters of a potential transit of HD 38529b. We 
present 12 years of precision photometry of HD 38529 ac- 
quired with an automated photoelectric telescope (APT) 
at Fairborn Observatory. These observations rule out a 
transit of the inner planet. Furthermore, we demonstrate 
that the host star exhibits a cyclic brightness variation 
on a timescale of approximately six years which is corre- 
lated with the S-index derived from our Keck I spectra. 

2. FUNDAMENTAL STELLAR PARAMETERS 

2.1. Stellar Radius 

HD 38529 was observed during three nights in Novem- 
ber 2012 using the Georgia State University Center for 
High Angular Resolution Astronomy (CHA RA) interfer- 
ometric array (ten Brummc laar et aI1l2005f) . Our obser- 
vational methods and s trate gy are described in det ail in 
Ivon Braun et al.l (|2012t ) and Boyaiian et~aT1 ()2012bO and 
references therein. We used two of CHARA's longest 
baselines (S1E1 and E1W1) to perform our observations 
in J?-band with the CHARA Classic beam comb iner 
(|Sturmann et al.l 120031 : Iten Brummelaar et al.l 120051) in 
single-baseline mode. We obtained 3, 1, and 2 observa- 
tions (brackets) during the nights of 3, 4, and 12 Novem- 
ber 2012, each of which contains approximately 2.5 min- 
utes of integration and 1.5 minutes of telescope slewing 
per object (target and calibrator). To remove the influ- 
ence of atmospheric and instrumental systematics, inter- 
ferometric observations consist of bracketed sequences of 
object and calibrator stars, chosen to be near-point-likc 
sources of similar brightness as HD 38529 and located at 
small angular distances from it. We originally used both 
HD 37077 and HD 36777 as calibrators but eliminated 
the latter due to the presence of a second fringe packet 
in each observation, indicating that HD 36777 may be an 
unresolved binary. 

To protect against unknown systematics in inter- 
ferometric data, we ordinarily require the use of at 
least two calibrators, two base lines, and data obtain ed 
during at least two nig hts dBovaiian et al.l l2012al lbl; 
Ivon Braun eTaD l2011aM 120121 1 Thus, we added 12 
archival CHARA AT-band brackets obtained in 2005 us- 
ing the S 1E1 baseline and H D 43318 as calibrator, pub- 
lished bv lBaines etaLl (|2008l) . 

The uniform disk and limb-darkened angular diameters 
($ud and #ld, respectively; see Table [J) are found by fit- 
ting the calibrated visibility measurements (Figure [lj to 
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Fig. 1. — Calibrated visibility observations along with the limb- 
darkened angular diameter fit for HD 38529 (top panel) along with 
the fractional residuals around the fit (bottom panel). The blue 
round points are the i-f-band data, and the red square s are the H- 
band data. The _ft"-band data from Baincs et al. (2008) were taken 
with a shorter effective baseline than our new if-band data (due 
to the difference in w avele ngth) . The reduced x 2 of this fit is 0.49. 
For more details, see i|2.1l 

the respective functions for each relation. These func- 
tions may be described as n th -order Bessel functions of 
the angular diameter of the star, the projected distance 
between the two telescopes, and the wavelength of ob- 
servat ion (see equations 2 and 4 of lHanburv Brown et al.l 
(|1974[ ) ) . Visibility is the normalized amplitude of the cor- 
relation of the light from two telescopes. It is a unitless 
number ranging from to 1, where implies no corre- 
lation and 1 implies perfect correlation. An unresolved 
source would have a perfect correlation of 1.0 indepen- 
dent of the distance between the telescopes (baseline). A 
resolved object will show a decrease in visibility with in- 
creasing baseline length. The shape of the visibility ver- 
sus baseline is a function of the topology of the observed 
object (the Fourier Transform of the object's shape). For 
a uniform disk this function is a Bessel function, and 
for this paper, we use a simple model for limb-darkening 
variation of a uniform disk. The visibility of any source is 
reduced by a non-perfect interferometer, and the point- 
like calibrators are needed to calibrate out the loss of 
coherence caused by instrumental effects. 

We use the linear limb-dark ening coeffici ent [Ih = 
0.362 from the ATLAS models in lClareti ([2000) for stellar 
T ff = 5500 K and log g — 4.0 to convert from £>tjd to 6*ld- 
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The uncertainties in the adopted limb darkening coeffi- 
cient amount to 0.2% when modifying the adopted grav- 
ity by 0.5 dex or the adopted T e g by 200K, well within 
the errors of our diameter estimate. 

Our interferometric measurements yield the following 
values for HD 38529's angular diameters: #ud = 0.593 ± 
0.016 milliarcseconds (mas) and 6*ld = 0.611 ±0.016 mas 
(Table [TJ. Combined wi t h the direct distance measure- 
ment from Ivan Leeuwenl (12007ft of 39.277 ± 0.617 pc, we 
derive a stellar radius for HD 38529 of 2.5780±0.0795i? Q , 
which is consistent with 2.44 ± 0.22i? Q calculated by 
iBaines et all (f2008h . 

2.2. Stellar Effective Temperature and Luminosity 

To calculate HD 38529's effective temperature and 
luminosity, we produce a spectral energy distribu- 
tion (SED) fit based on the spectral templates in 

M (11998f) t o pho to metry from Uohnson et al. | 
[Cousins! (fl962ft. lArgud (11966ft. iMermilliodl 
lHauck fc Mermilliodl (11998ft. lOlsenl (11993ft . 
Cutri et all (|2003t ). iMcClure fc Forrester! (|198lh . and 
Golavl (119721) . We furtherm ore use the distance calcu- 
lated in Ivan Leeuwenl (|2007l ) and set interstellar redden- 
ing to zero. 

From the SED fit, we calculate the value of HD38529's 
stellar bolometric flux to be F BO l = (12.02 ± 0.084) x 
10~ 8 erg cm~ 2 s _1 and, consequently, its luminosity L = 
5.777 ± 0.186Lq. The calculated effective temperature 
for HD 38529 is T cS = 5576 ± 74 K (see Tabled}. 

3. REFINING THE PLANETARY ORBITS 

Here we present the new RV data and the revised or- 
bital solution for the HD 38529 system. We combine 
this with the derived host star properties to determine 
an accurate transit ephemeris for HD 38529b. 

3.1. Spectra Acquisition 

The RV data for HD 38529 presented here comprise 
436 measurements and were acquired from 3 instru- 
ments/ tele scopes: the High Resolution Spectrograph 
(jTullI 119980 on the Hobby-Eberl y Telescop e (HET), the 
Hamilton Echelle Spectrograph (jVogtlll987l ) on the 3.0m 
Shane Telescope at L ick Observatory, and the HIRES 
echelle spectrometer (jVogt et al.l 11994ft on the 10.0m 
Keck I telescope. Shown in Table [2] are a subset of the 
full dataset, available in the electronic version of this 
paper. The fourth column in Table [5] indicates the six 
independent datasets that form the combined data, two 
each from HET, Lick, and Keck. The data sou r ces ar e 
as follows: 1 = HET data from Benedict et al.1 (|2010ft : 
2 = new HET data presented here; 3 = K eck data from 
Wrig ht et al.1 (|2009l ): 4 = Keck data from iWright et al.1 
(2001) before 2009 September 15, new K eck data pre- 
sented here thereafter; 5 = Lick data from IWright et 
( 2009); 6 = Lick data from IWright et al.l ((20091) before 
2009 September 15, new Lick data presented here there- 
after. The division of the Lick data into two separate 
datasets is necessitated by a change in the dewar result- 
ing in different CCD response characteristics. The offsets 
between these datasets are accounted for in the Keplerian 
orbital fitting described below. 

3.2. SME Analysis 



TABLE 2 
HD 38529 Radial Velocities^ 



Date RV a Datasett 
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13341.779899 


-105.27 


7.77 


1 


13341.898484 


-118.43 


7.25 


1 


13355.845730 


-102.05 


7.34 


1 


13357.859630 


-105.27 


7.48 


1 


13358.724097 


-87.82 


7.11 


1 


13359.729188 


-82.07 


8.70 


1 


13360.849520 


-65.85 


7.80 


1 


13365.817387 


1.45 


7.73 


1 


13367.812640 


-20.41 


9.48 


1 


15095.967984 


33.82 


6.08 


2 


15115.906428 


-45.13 


6.09 


2 


15141.862009 


34.65 


6.05 


2 


15142.942977 


2.05 


5.64 


2 


15175.755713 


-66.65 


5.90 


2 


15175.758488 


-68.20 


5.60 


2 


15176.860583 


-44.69 


5.07 


2 


15182.725710 


38.24 


4.68 


2 


15185.724484 


-15.63 


5.29 


2 


10418.959317 


75.60 


1.32 


3 


10545.771238 


12.67 


1.35 


3 


10787.014317 


-69.45 


1.40 


3 


10807.061991 


-70.45 


1.30 


3 


10837.758229 


-111.32 


1.50 


3 


10838.784387 


-113.34 


1.43 


3 


10861.729653 


-30.86 


1.42 


3 


10862.725174 


-41.86 


1.46 


3 


11073.058843 


-64.92 


1.32 


3 


13750.837303 


32.41 


1.18 


4 


13750.837905 


32.48 


1.32 


4 


13750.839132 


27.87 


1.02 


4 


13751.878079 


49.46 


0.99 


4 


13752.884016 


47.48 


0.92 


4 


13753.895081 


33.16 


0.98 


4 


13775.747373 


-5.19 


0.91 


4 


13776.880150 


10.45 


0.87 


4 


14336.130035 


261.65 


1.11 


1 


11101.015625 


-101.58 


6.57 


5 


11101.035156 


-111.09 


6.40 


5 


11102.014648 


-88.30 


6.17 


5 


11102.032227 


-87.47 


6.68 


5 


11131.912109 


-91.36 


5.99 


5 


11131.930664 


-79.77 


7.06 


5 


11132.928711 


-67.59 


6.02 


5 


11154.808594 


-140.06 


6.89 


5 


11173.906250 


-104.89 


10.69 


5 


12267.809570 


143.48 


3.88 


6 


12267.818359 


150.79 


3.86 


6 


12298.721680 


103.71 


4.99 


6 


12298.751953 


113.82 


5.33 


6 


12298.785156 


105.67 


5.06 


6 


12298.815430 


93.78 


6.53 


6 


12299.695312 


110.60 


4.16 


6 


12335.657227 


206.29 


4.74 


6 


12534.970703 


31.96 


4.15 


6 



t Shown here is a subset of the data. The complete 
dataset contains 436 measurements and is available 
electronically. 

t Datasets 1 and 2, 3 and 4, 5 and 6 are from 
telescopes HET, Keck, and Lick respectively. 

For additional insight into the properties of 
the host star, we used Spectroscopy Made Easy 
(jValenti fe Piskunovl 11996ft to fit high-resolution Keck 
spectra of HD 38529. The methodology of this 
technique, including application of the wavelength 
int ervals and line da t a, ar e described in more detail 
by iValenti fc Fischerl (|2005l ). We further constrained 
the surface gravi ty using Yonse i -Yale (Y 2 ) stellar 
structure models (iDemarque et al.l 120041 ) and revised 
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Hipparcos parallax es (Ivan Leeuwenl 12007ft with the 
iterative method of iValenti et al.l ( 2009ft . The resulting 
stellar parameters are listed in Table [1] along with 
the directly measured parameters from Ej2j The SME 
derived parameters are effective temperature, surface 
gravity, iron abundance, projected rotational velocity, 
mass, radius, and age. These properties are consistent 
with a slightly metal-rich, mid-G sub-giant. There 
are other literature sources that hav e analyzed the 
abund an ces in HP 38529 , namely: iGonzalez et al] 
( 2001ft. iZhao et al.l (l2002Tl. iBodaehee et al.l (12005 



cuvillo 




et all (120061). iGilli et all (120061). iTakeda et all 
Delgado Mena et al.l (12010ft. iBrugamyer et al.l 
iKang et al.l (|201lD . IPetigura fc Marcvl (j2011ft 



These catalogs contain [Fe/H] values measured with a 
variety of telescopes, techniques, and stellar models, and 
the determinations span the range of 0.28-0.46 dex with 
a mean of 0.38 dex and a median of 0.4 dex. These find- 
ings a re consistent with the results of IValenti fc Fischer! 
(2005]), where the average error for [Fe/H] is ±0.05 dex. 

The stellar radius is an essential parameter for esti- 
mating the depth and duration of a planetary transit. 
Although the spectroscopic and interferometric derived 
tcmcperatures are consistent with each other, the radii 
are not. This difference is a result of the SME luminos- 
ity, which uses bolometric corrections. These corrections 
can sometimes be unreliable for evolved stars such as 
HD 38529 and, in this case, is underestimated compared 
with the luminosity derived from the SED fit and paral- 
lax. Since a similar T c g is found from the two techniques, 
this results in a reduced radius estimate. As described in 
H3.51 we adopt the interferometric radius for the purposes 
of estimating the predicted transit properties. 

3.3. Keplerian Orbital Solution 

We fit a two-planet Keplerian orbital solution to the 
RV data using the partially li nearized, least-squares fit- 
ting procedure described in iWright fc Howardl (2009) 
and estimated parameter uncertainties using the BOOT- 
TRAN bootstrapping routines described in iWang et al.l 
(12012ft . Table [3] lists the resulting parameters for the 
two-planet orbital solution. The data and orbital solu- 
tions are shown separately for the b and c planets in 
Figure [21 along with the residual velocities with respect 
to the best 2-planet orbital solution. The improved mass 
estimate for planet c combi ned with the F i ne Gu idance 
Sensor (FGS) astrometry of lBenedict et al] (|2010ft yields 
a true mass es t imate of 16.76 ± 0.11 Mj. As noted by 
iBenedict et al.l (|2010f ) . this lies within the brown dwarf 
mass regime. 

The fit required five additional free parameters due to 
the offsets between the six independent datasets. The off- 
sets with respect to the observations acquired with HET 
dataset 2 (the new HET data presented in this paper, 
see 331} arc -44.9, -63.92, -73.01, -117.2, and -123.1 
ms _1 for datasets 1, 3, 4, 5, and 6 respectively. As shown 
in Table [3l the Xred an< ^ rms scatter of the residuals are 
relatively large. The major contributor to the calculated 
X^ed vame are the Keck data whose uncertainties are sig- 
nificantly smaller than those associated with the Lick and 
HET data. Specifically, the Keck rms scatter of 8.6 m s _1 
exceeds its internal error, which has a median value of 1.3 
ms _1 . This indicates that there is a stellar noise com- 
ponent to the overall noise level which is not accounted 
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Fig. 2. — Top and Middle Panels: Radial velocity signal (black 
dots) induced by HD 38529b and c, respectively, and the best-fit 
orbital solution (dashed line). Error bars shown are internal errors 
for each observation. The radial velocity signal for each planet was 
extracted by subtracting off the best-fit orbital velocities of the 
other planet from the total observed RVs. Bottom Panel: Residual 
velocities with respect to the best 2-planet orbital solution. The 
red dots are for Keck data (dataset 3 and 4), the blue triangles 
are for Lick data (dataset 5 and 6), and the green squares are for 
the HET data (dataset 1 and 2). The typical size of internal error 
bars for each telescope (± median internal errors) are plotted on 
the upper right of this panel. 
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TABLE 3 

Keplerian Orbital Model 



False Positive Probability: 5% 
10% 



Parameter 



Value 



HD 38529 b 
P (days) 

T c a (JD - 2,440,000) 
T p *> (JD - 2,440,000) 

e 

K (ms- 1 ) 
w (deg) 
Mpsini (Mj) 
a (AU) 

HD 38529 c 
P (days) 



T v b (JD - 2,440,000) 



K (ms- 1 ) 
w (deg) 
M p sini (Mj) 
a (AU) 

System Properties 

7 (ms ) 

Measurements and Model 

N ohe 

rms (ms -1 ) 

Y 2 



14.30978 ± 0.00033 
15815.633 ±0.063 
12281.19 ±0.15 
0.259 ± 0.016 
56.81 ± 1.01 
93.3 ±4.1 
0.8047 ± 0.0139 
0.1278 ±0.0006 



2133.54 ±3.31 
12264.49 ± 6.43 
0.3472 ± 0.0057 
170.54 ± 1.12 
20.08 ± 1.14 
12.51 ±0.08 
3.594 ±0.018 



96.34 ± 3.59 

436 
11.76 
11.67 



a Time of mid-transit. 

b Time of periastron passage. 

for in the fit. The causes of the stellar noise, including 
pu lsations and star spot activity, are discussed in detail 
bv lBenedict et afl (|2010D . Another hypothesis is that of 
an additional planet with an amplitude of K > 5 ms _1 . 
We explain in the following section that our data do not 
support detection of such an additional planet. 

3.4. A Third Planet? 



iBenedict et al.l (|2010l ) utilize their results to speculate 
on evidence for a third planet in the system. Thus, we 
also consider this possibility from our analysis since our 
RV data comprise a substantially larger dataset. As re- 
ported by Benedict et al. (2010), a coplanar orbital 
solution is only stable if the third planet has a period 
within the window of [33,445] days and an eccentricity 
of e < 0.3, or a period larger than their RV data base- 
line (> 10 years). For this reason, we focused our search 
for the third planet within the period window of [33, 445] 
days and constrained the eccentricity to be < 0.3. 

We first searched for strong periodic signals in the 
residuals of the two-planet Keplerian solution by fit- 
ting sinusoids to the residuals at different periods within 
[33,445] days (with 0.4 day step in period). The results 
are plotted as the solid line in Figure [3] We then esti- 
mated the false positive probability to see if any of the 
strong peaks are significant. We define the false posi- 
tive probability for a peak with a certain amplitude K' 
as the probability that a signal with amplitude > K' is 
generated by the residuals by chance. We generated 1000 
sets of simulated residuals by scrambling the true residu- 
als (and their associated errors, with replacements), and 
then searched for the peak with largest amplitude within 
the P = [33, 445] day window for each of the 1000 scram- 
bled data sets. These 1000 amplitudes provide approxi- 
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Fig. 3. — Amplitude of best-fit sinusoids to the residuals of the 
two-planet Keplerian solution (solid line). Any peak in this pe- 
riod window that has amplitude larger than the top dashed line is 
considered to be significant for having < 5% false positive probabil- 
ity. The two lower dashed lines (< 10% and < 50%) have similar 
meanings. No period within this window has less than 5% false 
positive probability, and the two peaks with < 10% false positive 
probability are at 119 days and 164 days. We see no sign ificant 
peak around 194 days as reported by Benedict ct al. (2010). See 
Section 3.4 for more details. 

mately the distribution of amplitudes arising purely from 
random noise in the residuals. Any peak in Figure[3]that 
has an amplitude smaller than 950 (95%) of these 1000 
amplitudes is thus considered having false positive prob- 
ability of > 5%. This is marked by the top dashed line 
in Figure [3] and similarly for the 10% and 50% lines. 

As shown in Figure [31 no peak has a false positive 
probability of less than 5%, and there are just two with 
less than 10% at 119 days and 164 days. The highest 
peak at 164 days has a false positive probability of 6.8%. 
W e see no significant pea k around 194 days as reported 
by IBenedict et al.l (|2010( ) . We then performed 3-planet 
Keplerian fit with our RV data within the P — [33, 445] 
day window and with the constraint that the eccentricity 
must be smaller than 0.3. We found that indeed the best- 
fit is near 164 days, with e = 0.3 (also true if we force the 
third planet to be on a circular orbit; best-fit e — 0.99 if 
no constraint on e is required). The xt of this fit is 9.58, 
and an f-test suggests that the 3-planet model provides 
a better fit though having 5 more parameters. However, 
the rms for this fit is 11.92 m s -1 , i.e., adding a third 
planet does not reduce the rms of the fit. Combining 
with the fact that this signal at P = 164 days does not 
have lower than 5% false positive probability, we cannot 
conclude that our data have detected a third planet in 
the HD 38529 system. 

We note here that including or excluding this third 
planet does not affect our transit exclusion analysis in 
the following sections, because the changes in the orbital 
parameters for both HD 38529b or c, after adding the 
third planet, are smaller than their error bars listed in 
Table [3 

3.5. Transit Ephemeris Refinement 

From the stellar and planetary properties listed in Ta- 
bles [T] and [31 we compute a refined transit ephemeris for 
HD 38529b. We use the directly measured radius of the 
host star determined from our CHARA observations. We 
approximate the radius of the planet as R p — l.ORj, 
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based upon th e mass -radius relationship described by 
iKane fc Gelinol (|2012t ). These properties lead to a pre- 
dicted transit duration of 0.33 days a nd a predicted tran- 
sit de pth of 0.17% or 0.0018 mag. IKane fc von Braunl 
(2008) show how the probability of a planetary transit 
is a strong function of both the eccentricity and the 
argument of periastron. The periastron argument of 
HD 38529 is particularly well-aligned for transit prob- 
ability enhancement since it is quite close to the optimal 
angle of 90° . Combined with the eccentricity, this results 
in a transit probability of 12.8%. However, if the planet 
were in a circular orbit, the duration would be increased 
to 0.43 days but the probability would be reduced to 
9.5%. 

The calculation of the transit mid-point shown in Table 
[3] was performed with a Monte-Carlo bootstrap, which 
propagates the uncertainty in this orbital parameter to 
the time of the transit. The resulting uncertainty in the 
transit mid-point is 0.063days = 91minutes. As such, 
the transit window is dominated by the predicted dura- 
tion rather than the mid-point uncertainty, making it a 
suita ble candidate for photometric follow-up (|Kane et al.l 
2009). Thus, we compute a revised transit ephemeris 
from the orbital fit, which we utilize in the following sec- 
tions. 

4. PHOTOMETRIC OBSERVATIONS 

We have acquired 1106 photometric observations of 
HD 38529 on 992 nights between 2000 November 28 
and 2012 March 31, all with the Til 0.80 m Auto- 
mated Photoelectric Telescope (APT) at Fairborn Ob- 
servatory in Arizona. The Til APT, one of several such 
telescopes operated at Fairborn by Tennessee State Uni- 
versity, is equipped with a two-channel precision pho- 
tometer that uses a dichroic filter and two EMI 9124QB 
bi-alkali photomultiplicr tubes to separate and simul- 
taneously measure the Stromgren b and y pass bands. 
We programmed the APT to make differential bright- 
ness measurements of our program star HD 38529 (P, 
V = 5.95, B - V = 0.77, G4 IV) with respect to the two 
comparison stars HD 38145 (CI, V = 7.89, B — V = 0.33, 
F0 V) and HD 40259 (C2, V = 7.86, B-V = 0.38, F0 V). 
From the raw counts in both pass bands, we compute the 
differential magnitudes P - CI, P - C2, and C2 - CI, 
correct them for atmospheric extinction, and transform 
them to the Stromgren system. To improve our pho- 
tometric precision, we combine the differential b and y 
observations into a single (6 + y)/2 passband, which we 
indicate with the subscript by. Furthermore, we compute 
the differential magnitudes of HD 38529 against the mean 
brightness of the two comparison stars. The resulting 
precision of the individual P — (CI + C2) /2t lJ differential 
magnitudes ranges between ~ 0.0010 mag and ~ 0.0015 
mag on good nights, as determined from the nightly scat- 
ter in the C2 — CI differential magnitudes of the constant 
comparison stars. Further details of our automatic tele- 
scopes, precision photometers, and o bserving and data 
reduction procedur es can be found in iHenryl ( 19991 ) and 
lEaton et all ((2003) and refer ences therein. 

The resulting 1106 differential magnitudes, spanning 
12 consecutive observing seasons, are summarized in Ta- 
ble S] and plotted in the top panel of Figure HJ after 
normalization to bring the seasonal means to a common 
level indicated by the horizonatal line in the top panel. 
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Fig. 4.— Top: The 1106 individual P- (CI + C2)/2 by differential 
magnitudes of HD 38529, acquired with the Til 0.8 m APT during 
12 observing seasons between 2000 and 2012. The observations in 
each observing season have been normalized to give each observ- 
ing season the same yearly mean. Middle: The 1106 observations 
phased with the orbital period and time of transit of companion 
b. The semi-amplitude of a least-squares sine fit to this orbital 
phase curve is 0.00014 ± 0.00007 mag, consistent with the absense 
of periodic light variability in HD 38529 on the radial velocity pe- 
riod. This provides strong confirmation of the existence of planet b. 
Bottom: The observations near phase 0.0 plotted on an expanded 
scale. The solid curve shows the estimated depth (~ 0.0018 mag) 
and duration (±0.012 phase units) of a central transit of compan- 
ion b. The ±lcr extent of the transit window is indicated by the 
vertical dotted lines. Transits are ruled out to a depth of ~ 0.0004 
mag. 

The normalization removes season-to-season variability 
in HD 38529 perhaps caused by a starspot cycle in this 
mildly active star (see below). The scatter in the nor- 
malized data from their grand mean is a = 0.0019 mag 
(standard deviation). This is slightly larger than our 
typical measurement precision given above, suggesting 
low- amplitude, night-to-night variability in HD 38529. 

The normalized observations from all 12 observing sea- 
sons are replotted in the middle panel of FigureHJ phased 
with the 14.3 day planetary orbital period of HD 38529b 
and the time of mid transit (T c ) from Table [3l A least- 
squares sinusoidal fit to the phased data gives a formal 
semi-amplitude of just 0.00014 ± 0.00007 mag, which 
limits any periodic brightness variability of the star on 
the orbital period to a very small fraction of one milli- 
magnitude (mmag). This rules out the possibility that 
the 14.3-day radial velocity variations are simply jitter 
induced by stellar activity, as has been d ocumented in 
slightl y more active st a rs, fo r insta nce, bv iQueloz et al.l 
([200lh . iPaulson et all (|200l . and iBoisse et al.l (|2012D . 
Instead, the lack of photometric variability confirms that 
the radial velocity variations in HD 38529 result from 
true planetary reflex motion. 
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TABLE 4 

Summary of Photometric Observations for HD 38529 



Observing 




Julian Date Range 


Sigma 


Prot 


Full Amplitude 


< P — CI > 


< P - C2 > 


< C2 - CI > 


Season 




(HJD - 2,400,000) 


(mag) 


fdavs) 


(mag) 


(mag) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


2000-01 


47 


51877-51997 


0.00181 






-1.8254 ±0.0002 


-1.8033 ±0.0002 


-0.0221 ±0.0001 


2001-02 


88 


52169-52368 


0.00190 


38.5 ±0.5 


0.0024 ± 0.0005 


-1.8288 ± 0.0002 


-1.8061 ±0.0002 


-0.0226 ± 0.0001 


2002-03 


76 


52532-52732 


0.00219 


36.3±0.2 a 


0.0042 ± 0.0007 


-1.8303 ±0.0002 


-1.8076 ±0.0002 


-0.0226 ±0.0001 


2003-04 


95 


52894-53094 


0.00166 


36.2 ±0.3 


0.0029 ± 0.0004 


-1.8288 ±0.0002 


-1.8056 ±0.0001 


-0.0232 ±0.0001 


2004-05 


103 


53258-53462 


0.00138 


35.7 ±0.3 


0.0022 ± 0.0003 


-1.8259 ±0.0001 


-1.8033 ±0.0001 


-0.0225 ± 0.0001 


2005-06 


101 


53627-53827 


0.00175 






-1.8240 ± 0.0001 


-1.8017 ±0.0001 


-0.0223 ± 0.0001 


2006-07 


118 


53996-54194 


0.00177 


38.0±0.2 a 


0.0026 ± 0.0004 


-1.8246 ±0.0001 


-1.8015 ±0.0001 


-0.0231 ±0.0001 


2007-08 


110 


54370-54556 


0.00187 


37.3 ± 0.2 a 


0.0021 ± 0.0005 


-1.8265 ± 0.0002 


-1.8037 ±0.0001 


-0.0227 ±0.0001 


2008-09 


104 


54734-54919 


0.00198 






-1.8273 ± 0.0002 


-1.8047 ±0.0002 


-0.0226 ± 0.0001 


2009-10 


86 


55092-55284 


0.00189 


36.5 ±0.3 


0.0034 ± 0.0006 


-1.8275 ±0.0002 


-1.8042 ±0.0001 


-0.0232 ±0.0001 


2010-11 


95 


55459-55650 


0.00184 






-1.8253 ± 0.0001 


-1.8022 ±0.0002 


-0.0230 ± 0.0001 


2011-12 


83 


55830-56018 


0.00157 


37.2 ±0.4 


0.0024 ± 0.0004 


-1.8252 ± 0.0002 


-1.8019 ±0.0001 


-0.0232 ±0.0001 



a Periodogram analysis gave half of the quoted period, implying the 
photometric periods and their errors in these cases to get P ro t- 

The photometric observations within ±0.06 phase 
units of mid-transit are plotted with an expanded scale in 
the lower panel of Figure 0J The solid curve shows pre- 
dicted transit phase (0.0), depth (~ 0.0018 mag), and 
duration (±0.012 phase units) of a central transit, all 
computed from the stellar radius in Table [T] and the or- 
bital elements of HD 38529b in Table El The vertical 
dotted lines give the ±lcr uncertainty in the timing of the 
transit window, based on the uncertainties in the stellar 
radius and the improved orbital elements from Tables[T]& 
[3j respectively. Our data set contains 1084 photometric 
observations that lie outside the predicted transit time 
(solid curve); these have a mean of —1.81462 ± 0.00005 
mag. There are 22 observations that fall in transit; these 
have a mean of -1.81474 ± 0.00041 mag. The difference 
is our "observed transit depth," -0.00012 ±0.00042 mag 
in the sense that the mean of the transit points is slightly 
brighter than the mean of the out-of-transit observations 
but are, nevertheless, consistent with zero to four deci- 
mal places. As the lower panel shows, we have sufficient 
data around predicted transit times to rule out transits 
to a depth of ~ 0.0004 mag. 

It has been suggested bv lAnglada-Escude et alJ ()2010f ) 
that planetary systems in 2:1 orbital resonance can be 
mistaken for a single planet in an eccentric orbit. If 
planet b has an orbital period that is half of the de- 
rived ~ 14.3 day period, then it would have an even 
higher transit probability and thus a higher chance of 
being detected in the photometric data. We performed a 
search through the 12-year photometric data set for pe- 
riodic transits around 7.155 days and around 14.3 days. 
No transits were detected at or near either period. Al- 
though a dynamical model of a 2:1 resonance would lead 
to aperiodic transit times, the complete phase coverage 
of our photometric observations over a range of periods 
around 14.3 days leads us to conclude that such transits 
are ruled out. 

As noted above, the scatter in the normalized data 
set is slightly larger than the expected measurement 
precision, suggesting the presence of small starspots on 
HD 38529. Starspots on the photospheres of solar- type 
stars allow the possibility of direct determination of stel- 
lar rotation periods from rotational modulation in the 
visibility of the spots and the consequent variability 



star had spots on both hemispheres at those epochs. We doubled the 

in the star's brightness (see, e.g.. iSimpson et ail [2010). 
Spots can also produce periodic radial velocity variations 
that can mimic the presence of a planetary companion. 
Therefore, we performed periodogram analyses for each 
of the 12 seasons of normalized photometry plotted in 
the top panel of Figure [5] and, indeed, found very low- 
amplitude periodic brightness variations in 8 of the 12 
observing seasons. Similar analyses of the twelve C2 — CI 
seasonal datasets evinced no significant periodicity in the 
comparison stars. Figure [5] shows a sample frequency 
spectrum and phase curve for Season 5. Complete re- 
sults are given in columns 5 & 6 of Table fJJ The sea- 
sonal photometric periods scatter about their mean with 
a standard deviation of ~ 1 day. The weighted mean of 
the 8 photometric periods is 37.0 ± 0.4 days, which we 
take to be our best determination of the star's rotation 
period. This period is based on far more observations 
th an the preliminary pe riods of 35.7 and 31.6 days given 
by iFischer et all (|2003) and Benedict et""aH (120101 ) . re- 
spectively. The seasonal peak-to-peak amplitudes in col- 
umn 6 range between 0.002 to 0.004 mag, indicating spot 
filling factors of only a few tenths of one percent. Both 
the low level of spottedness and the fact that the stellar 
rotation period is distinctly different from the 14-day ra- 
dial velocity period and its harmonics demonstrates that 
stellar activity (spots and plages) is not the source of the 
14-day radial velocity period. 

If the inclination of HD 38529's rotation axis is near 
90° , so that V sin i approximately equals the equato- 
rial rotation velocity, V eqi then the stellar radius and 
projected rotational velocity from Table [1] result in a 
rotation period of ~ 37.5 days, essentially identical to 
our observed value of P ro t — 37.0 days. This implies a 
very high inclination for the stellar rotation axis and in- 
creases the probability of transits of HD 38529b, assum- 
ing the stell ar equatorial a n d pla netary orbital planes 
are aligned. iMoutou et alJ (|2011f ) have shown, for in- 
stance, that the orbital planes of the hot jupiters HAT-P- 
8b, HAT-P-9b, HAT-P-16b, & HAT-P-23b are all closely 
alig ned with the s t ellar rotation axis. In the same pa- 
per, |Mouto^eE!aD (1201 ID examined 37 exoplanetary sys- 
tems that have accurately measured spin-orbit angles 
and found that spin-orbit misalignment occurs primar- 
ily for stars with T e ff > 6300 K (spectral class F6 
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Fig. 5. — Top: A frequency spectrum of the 2004-2005 (Season 5) 
photometric observations of HD 38529. The best frequency occurs 
at 0.02803 ± 0.00027 cycles per day. Bottom: The 103 Season 5 
observations phased with the corresponding best period of 35.7 
days. The phase curve shows coherent variability with a peak- 
to-peak amplitude of 0.002 mag, which we take to be rotational 
modulation of photospheric spots. Eight of the twelve observing 
seasons exhibit similar modulation (see Table 



and hotter). Thus, the orbital geometry, the high stel- 
lar inclination, and the likelihood of spin-orbit alignment 
are all favorable for transits of HD 38529b, making our 
non-detection particularly disappointing. 

Finally, we examine long-term variability in 
HD 38529's Ca II H and K indices and APT pho- 
tometry to look for evidence of magnetic cycles that 
might induce apparent radial velocity variations and 
so mimic the presence of a long-period planet. Both 
H and K emission and brightness variability are 
good proxies for st e llar magnetic activity (see, e.g., 
iBaliunas et~all I1995L 119981 : lLockwood et al.1 120071 and 
references therein). In the top panel of Figure [6j we 
plot seasonal means of the Mount Wilson S-index 
derived from our K eck I RV spectra as desc ri bed i n 
IWright et al.l (l200l and I saacson fc Fischerl (f2010h . 
Despite a couple of observing seasons without Keck 
H and K measurements, clear variability is seen on a 
timescale of several years. In panels 2 & 3, we plot 
the 12 seasonal mean P — CI and P — C2 differential 
magnitudes (Table HJ) without any normalization, as was 
applied to the P — (CI + C2)/2 differential magnitudes 
in Figure H These two lightcurves show that HD 38529 
varies in brightness by several mmag from year-to-year 
with respect to both comparison stars CI and C2. The 
bottom panel of Figure [6] plots the yearly mean C2 — CI 
comparison star differential magnitudes. The number 
in the lower left corner of each panel gives the range 
of the seasonal means in magnitudes; the number in 
the lower right corner gives the standard deviation of 
the seasonal means with respect to their grand mean, 
shown by the horizontal dotted line in each panel. The 
standard deviation of the C2 — CI observations is only 
0.00038 mag, indicating that both comparison stars have 
excellent long-term photometric stability. Therefore, the 
short- and long-term variability evident in the P — CI 
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Fig. 6. — Long-term variations in Mt. Wilson S-index (Top) and 
brightness (Panels 2 and 3) measured with Keck I and the Til 
APT. The cycle timescale is ~ 6 yrs in all three timeseries. Long- 
term differential magnitudes of the two comparison stars CI and 
C2 (Bottom) show excellent stability to ±0.00038 mag), demon- 
strating that the variability in the P — CI and P — C2 light curves is 
intrinsic to HD 38529. The direct correlation of S-index and bright- 
ness is typical for stars of solar age and older. The weak magnetic 
activity in HD 38529 cannot produce the large-amplitude radial 
velocity variations attributed to planet c. 

and P — C2 light curves must be intrinsic to HD 38529. 

We see cyclic variation in both the H and K and APT 
observations plotted in Figure [6] Simple least-squares, 
sine-fit periodogram analyses of the H and K, P — CI, 
and P-C2 means give periods of 6.29 ±0.37, 6.21 ±0.43, 
and 6.12 ± 0.47 years, respectively, all indentical within 
their uncertainties. Futhermore, the brightness of the 
star and the strength of the H and K emission vary di- 
rectly in phase with each other. The cycle timescale and 
amplitudes of the H and K and photometric observations 
are typical of similar long-term cycles seen in a large sam- 
ple of solar-type stars being measured with the APTs 
(see IHenrvi [19991: lLockwood etall [20071 IHall et al.1120011 
and are analogous to solar irradiance variations measured 
for o ver three decades with space-based radi ometers (see, 
e.g.. iWatson fc Hudsonlll981t lFrohlichll2012D . The direct 
correlation of long-term brightness and magnetic activ- 
ity in HD 38529 is in the same sense as the correlation 
between solar irradiance and S-index. Despite the simi- 
larities in age and spectral class of the Sun (4.6 Gyr, G2) 
and HD 38529 (4.4 Gyr, G5), the amplitudes of the S- 
index and brightness variability are larger in HD 38529 
by factors of approximate ly 2 and 5, respectively (see, 
e.g. lLockwood et al.l [20071 their Fig. 2). This is due to 
the fact that HD 38529 is a more massive, evolved sub- 
giant (Table 1) with a deeper convection zone. 

The radial velocity period of HD 38529c is 2134.9 days 
(Table [3]) or 5.85 years, very close to the spot-cycle 
timescale of 6.2 years. This raises the question that 
the radial velocity variations attributed to HD 38529c 
may, instead, originate from the stellar magnetic cycle. 
As mentioned above, starspots in active stars have been 
found to induce radial velocity variations on stellar rota- 
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tion timescales. Over the past 15 years or so, several au- 
thors have discussed the mechanisms by which long-term 
(decadal) magnetic cycles might be th e source of long- 
term radial velocity variations (e.g., Saar fe Donahue! 
19971 iSaar fe Fischerll2000USantos et al.ll2010f ). Recently. 
Lovis et aLl (|2011[ ) analyzed a sample of 304 FGK stars 
from the HARPS high-precision planet-search sample, 
looking for possible correlations of stellar activity cy- 
cles with radial velocity and spectral line-shape param- 
eters. They confirmed that stellar magnetic cycles can 
induce long-period, low-amplitude radial velocity varia- 
tions with amplitudes up to 25 ms _1 in the more active 
stars. For HD 38529c, however, its radial velocity semi- 
amplitude of 171 ms _1 makes it immune to this type of 
false positive. 

5. CONCLUSIONS 

With the discovery of extrasolar planets, characteriza- 
tion has become an important pursuit and is intrinsically 
linked to the understanding of the host star properties. 
Thus, even relatively bright stars, such as HD 38529, are 
worthy of further scrutiny in order to gain insight into 
properties of the planets hosted by the star. Here we have 
performed just such a task by presenting the most accu- 
rate properties of the star HD 38529 so far produced. 
Direct measurements of stellar radii via interfcromctry 
provide an essential test of stellar models and validation 
of radii derived from spectroscopy. In this case, we mea- 
sure a radius of 2.578 R© for HD 38529, which is con- 
sistent with a slightly metal-rich mid-G sub-giant. The 
new RV data for the star are used to calculate a new Ke- 
plerian orbital solution for the system that enables us to 



place significant limits on a previously postulated third 
planet. Combining our mass esti mate for the c compo - 
nent with the FGS astrometry of iBenedict et all (2010) 
allows us to confirm that this object does lie within the 
brown dwarf mass regime. The refined transit ephemeris 
from the Keplerian orbital solution is combined with 12 
years of precision photometry to demonstrate the vari- 
ability of the host star and the dispositive null detection 
of transits for HD 38529b. The importance of this study 
is clear since the combination and inter-dependence of all 
these effects leads to a greatly improved understanding 
of the system as a whole. The TERMS project is sys- 
tematically proceeding to provide such characterization 
of the planetary systems known to orbit the brightest 
stars in the sky. 
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